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摘要 研究 了 动态 围 堵 嫌 犯 问题 , 假设 网 络 边 长 相等 , 交 巡 警 与 嫌犯 的 速度 相等 . 建立 了 嫌犯 移动 
信息 更 新 下 的 交 巡 警 调度 问题 的 0-1 线 性 整数 规划 模型 , 模型 利用 点 截 集 条 件 使 调度 后 的 警力 形成 国 
堵 圈 , 并 对 嫌犯 的 逃跑 行为 建 模 , 由 此 得 到 了 动态 围 堵 嫌 犯 问题 的 动态 模拟 模型 . 算 例 考虑 分 割 非 等 
边 长 网 络 的 边 , 然后 将 分 割 后 的 网 络 视 为 等 边 长 网 络 . 


关键 词 ”网 络 优化 动态 围 堵 模 型 ， 交 巡警 服务 平台 0-1 整 数 规划 


1 引言 


文献 四 指出 2011 年 全 国 大 学 生 数 学 建 模 竞 赛 B 题 的 围 堵 嫌 犯 问题 仍 是 值得 进一步 研究 的 课题 , 动 
态 围 堵 嫌 犯 问题 也 是 值得 进 一 Rn 制 的 方法 研究 了 
网 格 图 形 上 多 警察 协作 围墙 系统 . 不 同 于 文献 四, 本 文 对 一 般 路 网 上 的 围 堵 嫌 犯 问 题 通过 数学 建 模 的 
方法 开展 初步 的 研究 . 

动态 围 堵 嫌犯 问题 中 交 巡 警 和 嫌犯 在 网 络 上 的 位 置 为 连续 变量 , 不 能 确保 同一 时 刻 交 巡警 和 嫌犯 
都 在 网 络 节 点 上 . 本 文 假设 网 络 边 长 相等 和 交 巡 警 与 嫌犯 的 速度 相等 , 从 而 可 以 用 离散 优化 的 方法 建 
模 . 建立 了 嫌犯 逃跑 信息 更 新 下 的 交 巡 警 重新 调度 的 0-1 线 性 整数 规划 模型 , 通过 对 嫌犯 的 逃跑 规则 的 
建 模 , 建立 了 动态 围 堵 问题 的 模拟 模型 . 在 算 例 中 , 对 文献 3] 和 4 的 算 例 的 网 络 的 边 进行 分 割 , 然后 视 
为 等 边 长 网 络 . 这 一 做 法 也 适合 解决 实际 问题 . 

文献 [四 给 出 了 点 截 集 判 断 的 优化 模型 , 并 用 该 模型 建立 了 围 堵 嫌 犯 问题 的 优化 模型 . 文献 [3] 和 [4 进 
一 步 考虑 了 交 巡 警 可 以 占据 某 些 路 口 节点 , 以 使 得 嫌犯 到 达 其 它 点 的 时 间 更 长 的 堵截 嫌犯 问题 . 本 文 
的 模型 利用 了 文献 [ 引 给 出 的 点 截 集 判 断 的 建 模 方法 , 以 使 调度 后 的 警力 形成 围 堵 圈 . 


2 ” 单 步 信息 下 的 模型 


2.1 ” 交 这 警 调 度 


无 向 网 络 G = (V,), 其 中 点 集 V = {1,2,… ,n}, 五 为 边 集 . 任意 边 (i, 7) E 的 边 长 dij = 1, 交 巡 
警 和 嫌犯 的 速度 都 为 1, 即 单位 时 间 内 交 巡 警 和 嫌犯 都 正好 移动 一 个 边 长 的 距离 . 
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实际 路 网 的 边 长 不 相等 , 可 以 将 边 长 为 di; 的 边 (i, 门 等 分 成 roud(qij/ 荆 段 , 其 中 工 为 常数 , roud 为 四 
售 五 入 函数 , L 越 小 ， 分 割 后 的 网 络 的 平均 边 长 越 小 ， 边 长 集合 的 方差 也 越 小 , 但 对 于 计算 来 说 问题 的 
规模 会 变 大 , 在 实际 问题 中 可 以 根据 需要 选择 的 大 小 . 将 分 割 后 的 网 络 视 为 等 边 长 网 络 . 

交 巡 警 初始 位 置 J c V, 记 交 巡警 平台 总 数 为 | 咱 . 一 个 交 巡 警 平 台 只 能 负责 一 个 路 口 . 在 围墙 过 程 
中 , 多 个 交 巡 和 警 平台 的 警力 可 能 和 同一 个 节点 相遇 

嫌犯 每 到 达 一 个 路 口 并 开始 朝 另 一 节点 行驶 的 时 刻 , 需 根 据 新 的 嫌犯 移动 信息 重新 调度 交 巡 和 警 以 
形成 更 好 的 围 堵 . 以 “0" 表 示 当 前 状态 “1 表示 指派 方案 执行 后 的 下 一 个 状态 . 已 知 嫌犯 的 当前 位 置 节 
点 so 和 嫌犯 的 行驶 方向 为 从 so 到 si, 边 (so, s1) < 五 , 称 嫌犯 移动 一 个 边 长 距离 的 信息 为 “ 单 步 信息 ”. 当 
前 节点 i 处 的 警力 为 99 个 交 巡 警 平 台 的 警力 . 需 确定 交 巡 警 指派 方案 p = {pij, (全力 € 五 }. 

指派 变量 pij 满 足 


pi; € [0, 0] NZ, (i,j) eb, 
其 中 Z 表 示 整 数 集 . 假设 pj; 和 pi; 中 至 少 有 一 个 为 0, 即 


pjipi; = 0, (i,j) EE. (1) 


Ql e [0, J NZ, ieV. 
警力 分 布 899、Q1 和 指派 方案 p 在 点 i 处 的 平衡 关系 表示 为 
Q@; 十 Spii— > Dj =QbieT (2) 
(j,iEBE (ij)EE 
假设 一 个 节点 只 需 一 个 交 巡 警 平台 的 警力 就 可 以 堵截 嫌犯 . 记 围 堵 变量 yk e {0,1}, 表示 在 警力 分 
布 Q* 下 节点 i 处 是 否 有 警力 , 1 表示 有 , 0 表示 无 . 警力 分 布 &9 为 已 知 量 , 围 堵 变 量 y 与 警力 分 布 & 的 关系 
表示 为 


Q* > ke{0,1},ieV; 
yl 2Q*, kef{f0,1},ieV. 


假设 嫌犯 当前 位 置 so、 行驶 目标 s1 和 当前 警力 分 布 8? 为 已 知 信息 , 需 依据 一 定 的 目标 和 约束 下 做 
出 交 巡 和 警 指派 决策 . 为 了 让 所 有 的 警力 朝 嫌犯 靠近 , 考虑 以 下 目标 
Tl(so, s1,Q0) = min >》， Qidis. 
iEV 
在 包围 圈 形成 后 , 该 目标 还 可 以 起 到 缩小 包围 圈 的 作用 . 
为 了 描述 交 近 警 对 嫌犯 的 围 堵 状 况 , 定义 0-1 变 量 ys 和 wi, 满足 


n+ <1l, ieV; 
i= ieV; 


其 中 wi; = 1 对 应 的 点 集合 表示 已 有 围墙 点 (y; = 1 对 应 的 点 ) 和 理想 的 新 增 围 堵 点 (yi = 1 对 应 的 点 ) 一 起 
构成 的 集合 ee i 二 1 对 应 的 点 集 将 点 s1 和 点 集 T 分 割 开 , 详 见 文献 [5]). 
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(1— wi)(l— wri= (1— wl -yw)r, (i E 五 ; 


>》 ri=0; 


ZE 人 


xi€E1{0,1}, iEV. 


以 包围 圈 缺 口 尽量 小 为 目标 可 表示 为 


Ta(s0,51,Q0) 二 min > ， wi. 


i€V 
当 目 标 值 rs = 0 时 , 警力 形成 有 效 围 堵 . 
综合 以 上 约束 得 约束 集 C1. 
约束 集 C1 
N+ 1 
pi = + 好， 
Tsl 1; 


(1 — wi) (1 — yi)ri = (1 — Yi)(l — wi)7;, 


>》 zi=0; 


iET 

QA) pi > 
(jieE (ij)EE 

Ql > 多 ， 

WIJ| > Qi, 

pjipi; = 0, 


pi; € [0,Q? NZ, 
Qi € [0, | NZ, 


Ti VY ,Yi wi € {0， 1}, 


iEV; 
1 E TY; 


(7I) EE; 


i€EV 


iEV; 

iEV 
(Gj) Ee Bi<j 
(7) € 五 ; 

1 ET; 

iEV. 


以 约束 集 C1 为 约束 和 ra 十 72 为 目标 , 得 到 交 巡 和 警 调度 模型 1 


模型 1 
T 一 Imin7rl 十 Tro2; 


二 :在 C1. 


引入 0-1 变 量 j4ij, 将 约束 (1) 线 性 化 为 


Di > Hij, (7) E 五 ; 

LijlJ| 2 pi (7) E 五; 
Ai thi 1 (人 四)E 五 < 

HUij E {0,1}, (7) E 成 


文献 可 给 出 了 (3) 的 线性 化 方法 , 因此 模型 1 可 化 为 0-1 线 性 整数 规划 模型 . 


2.2 ”模拟 流程 


假设 嫌犯 在 so 处 选择 下 一 个 点 s1 时 遵循 规则 1. 
规则 1: 记 嫌犯 的 路 径 中 so 的 前 一 个 点 为 s-1. 记 集合 


41 = {7|(s0, 7) € 已 ,QI 一 0,7 天 5_1}. 


若 41 关 ,嫌犯 在 该 集合 中 以 相等 的 概率 随机 选择 一 个 点 作为 1; 当 4 = 5 且 Q9 ，= 0 时 ,嫌犯 选 
择 。 ;作为 出 逃 点 , 即 嫌犯 折返 . 特别 的 , 当 嫌犯 和 交 过 警 同时 到 达 so, 嫌犯 仍 可 选择 s1, 继续 逃跑 
实际 问题 中 嫌犯 的 选择 是 不 可 控 的 , 实际 问题 中 可 以 由 专家 指导 或 根据 现 有 的 理论 分 析 嫌 犯 可 能 
的 行为 , 将 考虑 行为 的 目标 和 约束 加 入 模型 中 , 使 得 模型 更 有 应 用 价值 
假设 若 嫌犯 能 到 达 城 区 出 入 口 集合 C, 则 嫌犯 成 功 逃 跑 . 每 次 根据 嫌犯 的 信息 重新 指派 交 巡 警 , 模 
拟 计算 的 流程 为 流程 1. 


流程 1: 

wy 0 确定 初始 值 so, 8?. 初始 时 刻 不 存在 si 点 . 
5 1 若 41 = B 且 Q0 ,六 0, 嫌犯 被 捕 , 停止; 否则 转 下 一 步 

2 由 规则 1 确定 出 逃 点 si1. 若 sl € C, 嫌犯 逃脱 , 停止 ; 否则 转 
下 一 步 ; 
3 由 模型 1 求 最 优 指派 方案 下 的 @L. 

4 更 新 符号 s_1 := s0,s0 := s1),Q0 := Q1, 返回 1. 


3 多 步 信息 下 的 模型 


3.1 ， 交 巡警 调度 


不 同 于 前 一 节 , 本 节 假 设 原 网 络 为 G = (万 ), 分 割 后 的 新 网 络 为 G' = (V',B'). 记 新 网 络 的 边 e 在 
原 网 络 的 边 e 上 为 e CC e. 记 原 网 络 边 e e 上 的 所 有 的 新 网 络 上 的 边 构成 的 集合 为 B'(e) = {ele CC 
ee' € BP'}. 记 边 集 2'(e) 关 联 的 新 网 络 上 的 点 集 为 V'(B'(e)). 

嫌犯 当前 位 置 为 so, 选择 逃跑 目标 为 sj, 因 (so, s1) 被 分 割 为 多 条 边 , 称 此 时 嫌犯 从 so 逃 往 si 的 信息 
为 “多 步 信息 ”. 假设 嫌犯 的 决策 需 满足 边 (so, s1) 上 无 警力 . 为 了 简化 分 析 , 假设 嫌犯 选择 逃跑 目标 为 s1 
时 , 还 需 满足 没有 警力 可 以 先 于 嫌犯 到 达 s1. 这 种 假设 也 使 得 嫌犯 变 得 更 聪明 . 

多 步 信 息 更 新 下 两 次 更 新 时 间 间 隔 更 长 , 交 巡 和 警 移动 的 距离 可 以 更 长 , 因此 需 考虑 任意 点 ;处 的 交 
巡警 可 能 被 指派 到 其 它 任 意 点 j, 即 交 巡警 指派 方案 p = {pij,i,7 EV'}. 相对 于 模型 1, 约束 需 相 应 作 修 
改 . 例如 , 约束 (4-7) 修 改 为 


As 
Oo 
as 


pi biEV,; 


一 、 
CD 
ss 


Mithuia<1l, jEeV,i<) 
pia=0, iEV'; 
Hij E {0， 1}, 2 了 7 E VV. 
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约束 (2) 修 改 为 


Q@3 十 


JEV 


D2, Pi 2, P= Qh ieV, 


EW 


记 新 网 络 中 点 i 到 ;的 最 短路 长 度 为 qij. 原 网 络 和 新 网 络 的 相同 两 点 间 最 短路 长 度 相等 . 点 so 到 si 的 


其 中 jj; 与 指派 决策 pj; 的 关系 由 约束 (8) 
约束 集 C1 修 改 为 束 集 C2. 
约束 集 C2 


Hjidji < dsosi) 1,] € V,, 


和 (9) 给 出 , 0-1 变 量 yj; 表 示 是 否 


距离 dsos 已 知 , 等 于 边 (so, s1) 包 含 的 分 割 后 的 网 络 的 边 数 . 将 点 ;的 警力 指派 到 点 i 需 满 足 


节点 7 的 交 这 和 警 被 指派 到 节点 i 


(1 — wi)(1 


Wt <1 ET 
Wi=m+, iEV’,; 
Vr 起 
Wi)zi = 1 ph), eB; 
DE = 0; 
ieT 
+ > pi pi;=Q), iEeV,'; 
jeV’ jEV’ 
Qi 二 yi, iEV,; 
角川 > QT 
Pij2Hi 7 ET 
MWilJ|2 pi 7 ET 
Wijtui<1l, iIEV,i<) 
pi=0, iEV,; 
Miidji 所 dsos bj EV,; 
pi El0, 0] NZ, ijeV,’; 
Ql el[0,IJNZ, ieV’; 
Hi E {0,1}, ijeEV,; 
riyl,y, wie {0,1}, ieV. 


以 约束 集 C2 为 约束 和 mi 十 2 为 目 
模型 2 


标 , 得 到 交 巡 警 调度 模型 2. 


T 一 Imin7rl 十 Tro2; 


届 风 (2. 


3.2 ”模拟 流程 
嫌犯 在 so 选择 下 一 个 点 s1 时 , 需 满 足 边 (so, s1) 上 无 交 巡 警 , 即 
Q"(i) = 0; (10) 


iEV'(E'((s0,)))) 
需 满足 嫌犯 能 先 于 所 有 交 巡 警 到 达 节 点 7 九 即 

min{y dij,i € V'} > dso0j, (11) 
其 中 yw 表示 当前 状态 下 点 i 处 是 否 有 警力 . 由 于 假设 嫌犯 和 交 巡 警 的 速度 相等 , 条 件 (11) 满 足 时 (10) 一 
定 满足 . 知 假设 交 巡 警 和 嫌犯 的 速度 不 相等 , 特别 是 嫌犯 的 速度 更 快 时 , 则 条 件 (10) 和 (11) 需 同时 满足 . 
i 


A2 = {7|(s0,7) € E,min{y dij,i EV'} 2 ds07,7 # 8_1}. 

假设 嫌犯 在 so 处 选择 下 一 个 点 si 时 遵循 规则 2. 

规则 2: 记 嫌犯 的 路 径 中 so eV 的 前 一 个 点 为 s_1 € V. 各 42 产 @, 嫌犯 在 该 集合 中 以 相等 的 概率 随 
机 选择 一 个 点 作为 s1; 当 42 = B 且 min{y04qi;s_1,i EV'} > dsos_1 时 , 嫌犯 选择 s_; 作 为 出 逃 点 , 即 嫌犯 折 
返 . 特别 的 , 当 嫌 犯 和 交 巡 和 警 同 时 到 达 so, 嫌犯 仍 可 选择 sl, 继续 逃跑 . 
流程 2: 
0 确定 初始 值 s0, 80. 初始 时 刻 不 存在 s_: 点. 
1 若 45 = B 且 min{y0dis_,,i EV'} < ds。 ,嫌犯 被 捕 , 停止; 
否则 转 下 一 步 . 
2 由 规则 2 确定 出 逃 点 si. 若 sl1 e C, 嫌犯 逃脱 , 停止 ; 否则 转 
下 一 步 . 
由 模型 2 求 最 优 指 派 方案 下 的 Q@1. 
4 更 新 符号 s_1 := s0,s0 := 51, 8 := Q1, 返回 1. 


注 1: 若 知 道 嫌 犯 更 多 的 信息 , 例如 当前 位 置 为 so, 接 下 来 将 经 由 s', 然后 到 达 s1, (so s'), (s', s1) < 已 
这 种 信息 更 新 下 的 建 模 只 需 对 模型 2、 规 则 2 和 流程 2 稍 作 修改 . 


CD 


4 ”模拟 计算 


对 文献 和 肌 和 [5] 的 算 例 作 如 下 修改 , 将 边 长 为 dij 的 边 (i 站 等 分 成 roud(dij/D) 段 , 取 L = 2, 得 到 图 1 所 
示 网 络 , 其 边 长 集合 的 方差 为 0.0338. 在 计算 中 , 将 该 网 络 当 作 等 边 长 网 络 . 
由 模拟 流程 1 和 2 模拟 围 堵 嫌 犯 . 用 MATLAB+TYALMIP+Gurobi 求 解 模型 1 和 2. 读者 可 以 通过 电子 
邮件 联系 作者 获取 相关 程序 . 模型 1 的 指派 变量 个 数 为 2|E|, 而 模型 2 的 指派 变量 个 数 为 2|V|?, 单 步 信 
息 更 新 系统 模拟 的 每 一 步 迭 代 所 需 的 计算 时 间 更 短 . 通过 多 次 模拟 计算 , 发 现 单 步 信息 更 新 和 多 步 信 
息 更 新 系统 模拟 中 嫌犯 都 可 能 逃脱 , 也 可 能 被 捕 . 读者 还 可 以 考虑 其 它 目标 函数 , 并 通过 模拟 计算 寻 
找 更 好 的 目标 组 合 . 

单 步 信息 更 新 系统 模拟 中 , 目标 立 可 以 驱动 交 巡 警 朝 嫌犯 移动 , 但 目标 ro? 在 围 堵 的 初期 因为 一 步 
的 移动 距离 太 短 而 很 难 改善 目标 ra. 怎样 在 初期 让 交 巡 警 的 路 径 就 朝 形 成 有 效 围 圈 的 方向 发 展 , 从 而 
减少 嫌犯 逃脱 的 概率 , 是 值得 进一步 研究 的 问题 . 
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图 1 算 例 图 


Figure 1 Caption 


多 步 信 息 更 新 系统 模拟 中 , 嫌犯 选择 下 一 个 目标 点 时 以 原 网 络 为 准 , 交 近 警 指派 根据 嫌犯 在 原 网 
络 节点 的 逃跑 信息 更 新 重新 优化 , 对 于 交 巡 警 的 网 络 为 分 割 后 的 等 边 长 网 络 . 相对 于 单 步 信 息 更 新 系 
统 模拟 , 因为 原 网 络 的 边 都 为 分 割 后 的 多 条 边 组 成 , 每 次 更 新 指派 决策 时 警力 都 可 以 行驶 更 长 的 距离 ， 


目标 nz 可 以 更 好 的 发 挥 作用 . 而 假设 嫌犯 也 知道 交 巡 警 是否 能 先 于 自己 到 达 某 个 路 口 


了 解 了 更 多 的 交 巡 警 的 信息 , 从 而 更 容易 逃脱 . 


5 小结 


节点 , 使 得 嫌犯 


本 文 从 等 边 长 网 络 和 交 巡 警 速 度 相等 的 假设 开始 考虑 , 有 效 的 解决 了 动态 围 堵 嫌犯 问题 中 交 巡 和 警 


和 嫌犯 在 网 络 上 的 位 置 为 连续 变量 这 一 难以 建 模 的 问题 . 对 于 交 巡 警 和 硬 


题 是 值得 进一步 研究 的 课题 . 


kh 犯 的 速度 不 同 的 动态 围 堵 问 


模拟 计算 中 , 嫌犯 被 围 堵 在 很 小 的 范围 内 , 读者 也 可 以 考虑 将 形成 有 效 的 包围 圈 ( 即 x2 = 0) 当 作 围 


堵 最 终 目 标 . 读者 还 可 以 根据 进一步 对 实际 问题 的 研究 , 或 通 
型 , 以 使 得 模型 更 有 应 用 价值 . 


过 和 专家 研讨 , 寻找 更 好 的 目标 或 修改 模 


0 


前 关于 动态 围 堵 嫌 犯 问题 以 及 其 它 动态 网 络 优化 问题 
优化 及 其 应 用 的 研究 有 一 定 的 启发 作用 . 


的 下 


完成 果 很 少 , 本 文 的 有 


究 对 动态 网 络 
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Dynamic suspect encirclement model over network with equal 
length edges 


ZHOU WeiGang*, FENG QianQian & CHEN ShiJun 
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Abstract This paper studies traffic and patrol polices’ suspect dynamic encirclement problem. We assume that 
the network edges have equal length and the polices and suspect have equal speed. First, we develop a 0-1 linear 
integer programming model for the police assignment problem with suspect movement information update, which 
uses the vertex cut conditions to let the polices to form an encirclement. Then, we model the rule of the suspect’s 
escaping behavior. Last, we obtain a simulation model for the dynamic encirclement problem. We cut the edges 
into small parts, and then regard the new network as with equal length edges in a numerical example. 


Keywords network optimization, dynamic encirclement model, patrol service platform, 0-1 integer program- 
ming 


